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| We present a predictive scheme for fermion masses and mixings based on supersymmetric 

■ SO (10) which provides an understanding of the observed hierarchies in the charged fermion 

sector including CKM mixings. In the neutrino sector bi-large mixing can be realized, while 
i-G ' the third leptonic mixing angle #13 ~ 0.01 — 0.2. 

CD . 

The SO (10) gauge group has many attractive features for building Grand Unified Theories 
(GUT) [1]. The standard model (SM) chiral fermions together with the SM singlet right handed neu- 
^ . trinos (RHN) can be unified in the 16 dimensional spinorial representations: 16 D (q, u c , d c , I, e c , z/ c ). 
The RHN states help generate non zero neutrino masses which can account for atmospheric [2] and 
solar [3] neutrino oscillations. The masses of v c states arise from the breaking of lepton number 
and this ingredient proves decisive for generating the observed baryon asymmetry of the Universe 
through leptogenesis [4]. The supersymmetric (SUSY) version of 5*0(10) strongly suggests unifica- 
tion of the SM gauge couplings at scale M G ~ 2 • 10 16 GeV. 

Many attempts have been made for building phenomenologically realistic models based on 
50(10) [5]- [17]. The problems associated with 5O(10) model building arise both from the scalar 
as well as the fermion sector. The scalar sector should be arranged to yield a suitable symmetry 
breaking pattern with preferably (for MSSM unification) a one step symmetry reduction, 5O(10) — ► 
SU(3)c x SU (2) L x U(1)y = G 32 i- At the same time, it is challenging to have natural doublet-triplet 
(DT) splitting in order to avoid rapid nucleon decay, and also maintain successful unification of the 
three SM gauge couplings. For resolving the DT splitting problem the mechanism in [6] (often called 
the missing VEV mechanism) can be invoked. Although its realization requires a complicated higgs 
sector supplemented by additional symmetries needed to avoid unwanted operators [8], [12], [15], 
it is a nice feature of 5O(10) that DT splitting can be realized without fine tuning. 
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In this paper we do not consider the details of the scalar sector and concentrate mainly on the 
fermion sector. Therefore, our choice of the SO (10) non-trivial scalar states will be as minimal as 
possible. For breaking SUSY SO (10) down to G 32 i, it is enough to introduce the scalar superfields 
£(45) + C(16) + (7(16), where the brackets display the dimensions of the SO (10) representations. 
The rank reduction from five to four is achieved by the VEVs (C) = (C) in v c direction, while (S) 
can be in a suitable direction preserving G321. To have a renormalizable Yukawa coupling for the 
top quark, we also introduce H(10) state in which the light h u + hd pair of MSSM higgs doublet 
superfields reside. Note that, in principle, nothing prevents the C, C states from also contributing 
with some weights to the physical doublet states [7]. In fact, as we will see below, this can play an 
important role for realizing realistic pattern for fermion masses and mixings. 

Turning to the fermion sector, let us note that the minimal version of SUSY 5*0(10) (containing 
only three matter 16-plets and a single higgs 10-plet) gives the incorrect asymptotic relations 
M v oc M D = M E , V C km = 1 and is therefore excluded. To build a realistic model some extension 
is needed and numerous attempts have been considered (for an incomplete list see [5]- [17]). In this 
paper we wish to build a simple and economical scenario with realistic fermion pattern, and which 
also accounts for bi-large atmospheric and solar neutrino mixings. The task is not trivial and to 
achieve our goal we first formulate the principles which serve us as a guide for model building: 

1) The field content must be as economical as possible in order to keep the model simple and 
maintain perturbativity up to the cut off scale. 

2) We wish to have the cut off scale at or near Mp\ ~ 2.4 • 10 18 GeV, so there is a large enough 
gap between M G and the cut off. Therefore there is a large interval where 50(10) is unbroken, and 
with a fixed field content one can check perturbativity of the unified gauge coupling 0:50(10) up to 
Mpi, taking into account all relevant threshold corrections. 

3) All operators allowed by the symmetries must be taken into account, so that the conclusions 
are robust. 

4) Last but not least, we would like to make the selection of the states and symmetries in such 
a way as to limit the allowed couplings and maximize the predictions. The third leptonic mixing 
angle #13 is a particularly important target in this regard. 

In the fermion sector, in addition to the three 16-plets, one can also introduce fermionic 10-plet 
states. The latter contain vector-like states with the quantum numbers of the left handed doublet 
and down type quark. Because of this, the 10-plets can play an important role [12], [13], [14], in 
the generation of charged fermion masses and mixings. In addition, one may need to introduce 
(scalar or fermionic) singlet superfields. They can play an important role in generating the desirable 
hierarchies and in explaining neutrino oscillations. In order to protect the hierarchies and create 
predictive power, we introduce an i?-symmetry which will play a central role in our considerations. 

The scalar and fermionic states (involved in the charged fermion sector) are given in table 1, 
which also lists the i?-charges of the corresponding state. S and 5 denote scalar superfields which 
are singlets under 50(10), while 10 and 10' are matter supermultiplets. The role of each state and 
i?-charge assignment should become clear shortly. Under R, the superfield fa and the superpotential 
respectively transform as fa — > e 1Qi 0j and W — > e^ 2a3+aH ^W . Table 1 lists all non trivial 5O(10) 
representations, so we can immediately check whether perturbativity works or not up to Mp\. 

Assuming for simplicity that below M G we just have the MSSM field content, the additional 
states enter into play above M G) so that a~ 1 ( / u) = a~ 1 (M G ) + ^ In (for /i > M G ). Taking 
a(M G ) ~ 1/24.4 we get a(M P \) ~ 1/26.7. Therefore, we remain in perfect perturbative regime 
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Table 1: i?-charges of introduced states. 





£(45) 


C(16) 


(7(16) 


5 


5 


R 


a s 


|a:/r-|a;s+2a; s 


-|aif-|a5+2as 


a^ 






#(10) 


16i 


16 2 , 16 3 


10 


10' 


R 


a^ 


a 3 -3a s +|a5 


«3-«e, a 3 


a 3 +!a//-fas-a s 


a 3 +|aH+fa s -a s 



which ensures that the prediction for gauge coupling unification can be trusted and also all possible 
GUT threshold corrections are calculable. This allows one to introduce, if needed, additional 
superheavy vector-like states. For instance, up to 8 pairs of 16 + 16 would be allowed . To fix the 
allowed field content from the requirement of perturbativity up to M P1 , let us point out that it is 
possible to introduce several 45 and 54 states with masses > Mq- For example, the allowed numbers 
which preserve perturbativity are (n^, n^) = (5,0), (4, 1), (1,3). Several scalar 45-plets together 
with 54-plet can play an important role in the solution of DT splitting problem [8], [12], [15]. Larger 
representations such 126 and 210 (which can play crucial roles in symmetry breaking and fermion 
mass matrices [17]) are not favored from this viewpoint because the 5*0(10) gauge coupling blows 
up not far above the scale Mq, thus not allowing a sufficiently high cut off scale. We will not 
pursue a detailed study of these issues here but only mention that the selection of a high cut off 
scale M* pa M P i is also important for suppressing operators such as F^ u , which introduce 

unknown threshold corrections [18]. With M* pa M P i, these corrections can be safely ignored. 

The hierarchies between charged fermion Yukawa couplings and the CKM mixing elements can 
be parameterized by the parameter A ~ 0.2 as follows: 

A* ~ 1 , A u : A c : A t ~ A 8 : A 4 : 1 . 

A 6 ~ A r ~ — tan P , \d : \ s : A& ~ A 4 : A 2 : 1 . 

m t 

A e : A^ : A T ~ A 5 : A 2 : 1 . 

V us pa\, V cb pa\\ V ub = A 4 - A 3 , (1) 

where the MSSM parameter tan/5 = Our aim is to gain a natural understanding of this 
hierarchical pattern. As expansion parameters, in addition to A ~(0.2), we exploit the dimensionless 
parameter ~ A 2 ~ 10~ 2 . Powers of this naturally appear from the Planck scale suppressed non 
renormalizable operators. Note that we take (C) = (C) ~ (E) ~ M G . In addition we assume that 
the scalar components of the superfields S, S also develop VEVs close to M G . Thus, 

(S) (S) (C) (C) (E) _ 
Mpi Mpi Mpi Mpi Mpi 
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With the i?-charge assignments in table 1, the couplings of 16-plets with H(10) are 



16i / 

16 2 
16 3 



I61 

S 2 CC 
"Mp7 



V 








16 2 


I63 





\ 








s 


Mpi i 


Mpi 


E 
Mpi 


1 / 



(3) 



where in each entry a dimensionless couplings of order unity is understood. Due to the appearance 
of operators of varying dimensions, there will be hierarchical structure between masses and mixings 
of different generations. The selection given in table 1 reduces the number of parameters and allows 
one to make several predictions. 

From (3) the asymptotic up type quark Yukawas are 



A* 



X u : A c 



: At ~ e 



(4) 



Note that the contribution of SO(10) violating VEV (E) is strongly suppressed for the third gen- 
eration 163-plet, so that at Mq we have 



A* — Ah — A T 



(5) 



On the other hand, the effect of (E) in the couplings involving second and first generations is strong. 
This means that the rotations arising from up and down left handed quark flavors are misaligned. 
This gives non trivial contribution to the CKM mixing matrix. From (3) we see that V c b ~ e has 
the correct magnitude. However, V us and V u b vanish at this stage and the contributions to the first 
and second generation masses of charged leptons and down quarks are negligible. Thus, additional 
contributions are necessary. We now come to the important role played by 10, lO'-plets and C, C 
states 3 . The relevant couplings (allowed by SO(10) x R symmetry) are 



Mi 



■10- 10' + 



9 CC C 

' 10C16 2 + — - 16! • 16 2 + — - 107n6 2 



pi 



M 



pi 



M 



pi 



M 



(6) 



pi 



These four couplings completely determine the masses of light generation of down quarks and 
charged leptons as well as CKM mixing angles V us and V u j,. In order to study the effects of these 
couplings, it is useful to work in terms of 577(5) multiplets. With 16« = 10j + 5j + lj, 10 = 5 + 5, 
10' = 5' + 5', and substituting appropriate VEVs, the couplings (3), (6) yield the following mass 
matrix relevant for down quark and charged lepton masses: 







_5i 


5 2 


5 3 


5 


5' 




10i 


/ 


e 4 h H 


echc 











\ 


10 2 




£ch c 


e 2 h H 


eh H 


eh c 


eh H 




10 3 







eh H 


h H 










5' 













6(E) 







5 


V 





6(C) 








6(E) 


/ 



(7) 



3 One can check that an extension with only 10-plets or C, C states cannot yield a satisfactory picture. It seems 
necessary to invoke both 10, 10' and C, C states. 
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where hn and he denote SU(5) 5-plets coming from H and C superfields respectively. With 
e c (~ 10~ 2 ) we denote (C)/M P1 in order to indicate (where appropriate) that there is no SU(5) 
breaking effect in the corresponding coupling. We see that there is strong mixing between 5 2 and 
5' states and it involves the SU(5) violating VEV (E). This will ensure that the wrong asymptotic 
relations m s = m M , m e = m d (common for minimal 5*0(10) and 577(5) GUTs) are avoided. Here 
we also assume that the light down higgs doublet h d resides in H and C (in 5 H and 5 C in terms of 
SU(5)): 

H2h d , CD\h d . (8) 

Upon integration of 5, 5, 5' and 5' states, and keeping in mind the content of SU(5) plets and 
composition of their tensor products 10 D q, e c , 5 D d c , I, 10 x 5 x 5h,c = qd c h d + e c lh d + • • • , (7) 
yields the following 3x3 mass matrices for down quarks and charged leptons, 




d c 2 (k) 
e c X 

e 

e 




(9) 



From (6) we see that 16 3 does not couple with 10, 10' states and therefore the relations in (5) are 
not affected. However, we have now generated masses for the second generation states with the 
required hierarchies 

— (m,^m„). (10) 

nib m r 

The CKM mixing angles are 

V us ~ A , V uh ~ V cb X , with V cb ~ e . (11) 

The (i-quark and the electron acquire masses through mixings with their nearest neighbors 
respectively (see (9)). This yields the following relations 

m d = \V us \ 2 m s , m e = \V e ^\ 2 m^ . (12) 

There is no SU(5) symmetry violating effects in the (1,2) and (2, 1) entries of (9). This leads to 
another predictive relation at the high scale, 

— = — . (13) 
m d 

With ^ ~ 3 near M G [19], we have the desirable asymptotic relation ^ ps 0.3. Thus, with a 
relatively simple 5*0(10) model we have obtained a realistic pattern of charged fermion masses and 
their mixings. 

Finally, the leptonic mixing angles acquire small contributions from the charged lepton sector, 
to wit, 9 efJl ~ y/^, 9^ T ~ e. Therefore, the neutral lepton sector (involving right handed neutri- 
nos) must play an important role in order to obtain a satisfactory explanation of the solar and 
atmospheric neutrino anomalies. As we will see now the 10, 10' states introduced above are also 
important for recovering bi-large neutrino mixings. 
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The Dirac mass matrix for v, v c differs from the case when there are only fermionic 16j-plets. 
The light lepton doublets by appropriate weights reside also in the 10-plets. This is indeed crucial 
point for natural bi-large neutrino mixings in our model. The couplings in eqs. (3), (6) yield the 
following matrix 



(14) 



where hu denotes the up-type higgs doublet coming from H. Integration of k, ly states lead to 
the 3x3 mass matrix 







v\ 






h' 


/5 


h 


f 


e 4 h H 











\ 


k 







e 2 h H 


eh H 





6(C) 


h 







eh H 


h H 








h 













6(E) 





h' 


V 





eh H 











m D 




(15) 



We see that (15) differs in structure from the up type quark mass matrix. In particular, the (2,2) 
and (3, 2) elements in (15) are of the same order (in e). This offers a natural way for large — v T 



mixing. However, we should also ensure the generation of the required value of Am^ 



/Am 



atm 



1/5, 



i.e. at leading order only one eigenstate should have mass close to a/ Am^ tm 4 . This is possible 
to realize if effectively in the see-saw mechanism only one RHN exchange dominates (the single 
RHN dominance mechanism) [20], [21]. For this to be realized the other two RHN states should 
decouple [21] without significant contribution to the light neutrino mass matrix. 

To achieve this, we introduce three 5*0(10) singlet states iVj (i = 1,2,3) with i?-charges 

3 1 1 

aN 2 = «3 + -®h — «s — -a>s , ctiVi = ®n 2 + 2«s + -as , ®n s = un 2 ~ a s + «s • (16) 



Note that for the time being not all phases in table 1 and (16) are fixed. Taking 



1 

®h = «s + -a s , 



(17) 



the couplings involving iV-states are 






N 2 


N 3 








N 2 


N 3 


1 





















s 


s 

Mpi 




iV 2 





1 


I) 








1 J 




N 3 


yo 








Mi 



pi 



(18) 



From eqs. (15), (18) the 9x9 mass matrix for neutral fermions is given by 



V 


v c 


N 





m D 


A 


T 





M 


A T 


M T 


M N 



(19) 



4 This requirement must be fulfilled because the neutrino masses turn out to be hierarchical. 
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where M and A emerge from the first matrix coupling of (18), while Mjv accounts for iViV-type 
coupling. The A contribution is non zero if the up type higgs doublet kg from C has non zero 
VEV. Integration of the heavy z/ c , N states induces contributions to the light neutrino mass matrix 
given by the double see-saw formula 

m?) = m D ^M N ^m T D - A±m T D - m D ^A T . (20) 

To analyze this mass matrix it is convenient to work in a basis in which the first coupling matrix 
of (18) is diagonal. With this change of basis the structures of mu and Mjv will be unchanged and 
we can parameterize the relevant matrices as 

/ e 4 \ 
m D = ( (3e ae \ h H , 
\ ae 1 J 

= Diag (e, m, e) • — ^ , M N = Diag (0, n 2 , 0) M P1 , (21) 

where the dimensionless couplings a,j3,rii t 2 of order unity have been identified in order to demon- 
strate that our goal (for which iV-states have been introduced) is achieved. 

Employing (20), the dominant contribution to the 3x3 neutrino mass matrix is given by 



mi 0) = ( (3 2 a(3\ , with to = (\a\ 2 + \(5\ 2 )n\n 2 ^ . ( 22 ) 

We see that only one eigenstate of (22) is massive with m Vz = m ~ m"^ 2 = ®^ ~ ®^ e ^ " ^ ne 
scale relevant for atmospheric neutrino oscillations. Moreover, the 2-3 mixing angle 823 is naturally 
large for The form in (22) emerges from the first term in (20). The second and third 

terms of (20) induce (negligible) corrections irrelevant for neutrino oscillations. 

For explaining the solar neutrino anomaly, we introduce an additional singlet state M with the 
i?-charge ajy = a 3 + \olh + 2a s — \ols- Note that we still have some freedom in selection of phases 
(determining the i?-charges), and we use this freedom to take a^, = ^as,a-g = -p^cts- With this 
and the prescriptions in table 1, and noting (17), the relevant superpotential couplings are 

Smo'CC + S 2 NlQiC + S 7 AT 2 , (23) 

where we have omitted appropriate powers of Mpi which can be restored when needed. Note that 
the lO'-plet contains the 1% (1/2) state with weight ~ 1. So, the coupling of 10' in (23) is of great 
importance. With this setting the coupling STiAfWsC is also allowed. Since this coupling is relevant 
we will take it into account for generality. Substituting appropriate VEVs and integrating the Af 
state we get a sub dominant contribution to the light neutrino mass matrix 




TO, 



otj3 ccy 

« = ( aP P 2 M ) to' , ( 24 ) 




Pi 7 



7 



where the dimensionless couplings a, f3, 7 ~ 1 and 

- ■ (25) 

We already see that the mixing angle 9 12 is naturally large for a ~ /3 ~ 1. As far as the mass scale 
m' is concerned, assuming the light physical up type higgs doublet is contained in H and C with 
weights 

HDK, CD sJUi u , (26) 

we have m! — 2 • 10 -3 — 5 • 1CT 2 eV, which accounts for solar neutrino oscillations. 

Next we study in detail the full neutrino mass matrix. Although in the neutrino sector there 
were many unknown parameters, the effective 3x3 neutrino mass matrix has a relatively simple 
structure. Because of this, the model enables one to calculate the third mixing angle 9\% in terms 
of other measured oscillation parameters. Once more we emphasize that this is possible thanks to 
the fermionic 10-plets. 

The neutrino mass matrix can be written as 

m v = mf + m« , (27) 

where denotes the dominant part responsible for the atmospheric anomaly, while the sub- 
dominant entry ensures large angle solar neutrino oscillations. They are given in (22) and 
(24), with m'/m <C 1. The leading part can be diagonalized by the transformation £/|jm^[/23 = 
m< 0) = Diag(0, 0, m), with 

/ 1 , , 

U 23 =Pl0, C23, S23 ) , ( 28 ) 

\0 , -s 2 3 , c 23 
where c 2 3 = cos #23, s 2 3 = sin #23, and 

tan ^3 = |4 , P = Diag (l, e 1 *, l) , X = -Arg(a/3) . (29) 
\a\ 

By this rotation, the sub-leading part transforms into 

(a 2 af3 0:7 \ 
a(3 (3 2 fa • w! , 
«7 Pi I 2 J 

with = (^c 23 e ix - 7s 23 ) , 7 = (/3s 23 e ix + 7C23) • (30) 
From (30) we see that for a ~ (3 the angle 9 12 is naturally large and therefore bi-large neutrino 

/ /~~A — ^ — 

mixing is realized. On the other hand, the third mixing angle 9 is ~ — ~ \ A m !° l is properly 

in y ZAm atm 

suppressed 5 . 

5 A similar relation was obtained in a democratic scenario [22] with discrete symmetries. We note, however, that 
so far it has not been generated within a GUT scenario supplemented by symmetry arguments. 
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The model allows one to make one more prediction if somehow the coupling Afl6 3 C can be 
forbidden. This is realized with a suitable modification of P-charges. For example, if the last 
coupling in (6) is generated through the integration of some additional states. Introduce the 
vector-like fermionic pair P(16) + P(16) and, instead of the last term in (6), consider the couplings 

16 2 PP + -f-10'FC' + SFF . (31) 
Mpi 

One can verify that with substitution of appropriate VEVs and upon integration of F, F states, 
the operator H1Q 2 ~ j^lO'-^l^ ls g enera ted, so that all results obtained in the charged 

fermion sector are robust. However, with the couplings in (31) the corresponding P-charges are 
modified as follows: 

13 1 

«c «c - ^("s - otg) , a v ^ a v - -(a s - otg) , a i6l -> a i6l + -(as - a-g) , 

113 

aio -> aio + -(a 5 - ag) , aio' — > aio' -|(«s- ag) , a Ni -> a^ + -(as - ag) , 

a a/" *■ a^ + -(as - ag) . (32) 

The P-charges of states which are not given here are unchanged. With this changes and with the 
following relations between phases: 

1 3 79 3 

otn = a s - -as + -oig , a s = — a s , agr = --ag , (33) 

all of the couplings presented above, with exception of jVl6 3 (7, survive. With 7 = 0, (30) reduces 
to 



m 



a 2 _ c 23 aj3 s 23 af3_ 
(i) = U^U 23 = P • c 23 c^ c 2 # c 23 s 23 /3 2 • Pm' , ( 34 ) 



«23a/3 c 23 s 23 /3 2 s 2 ,^ 2 



where P = Diag(l, e _1 ' x , e _1 ' x ). 

The matrix m u = mffi Imj, 1 ' is diagonalized by the transformation U^ 2 U^m v Ux 3 U 12 = m^ iag , 
where 

/ cia , sia , \ /c 13 e i5 , , s^e 1 - 5 

C/ 12 ~ P' • -sia , c 12 , , U 13 ~ , 1 , 

V , 0,1/ V - fi 13 , , C13 

P' = Diag (e^', e" 1 * e" 1 *) , X ' = -Arg(a0) , 5 = -Arg (j^j , (35) 

s 23 |a/5| , (36) 



m' 



tan 6>i2 = — , tan 6>i 3 
m^ iag = Diag(0,m 2 ,m) , m 2 = m'(|a| 2 + c 2 3 |^| 2 ) . (37) 
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For the 0i 3 mixing angle we have introduced the superscript V in order to indicate that this 
contribution comes from the neutrino sector. As we will see, the whole leptonic mixing angle 
9 13 = \U l e3 \ receives sizable contribution also from the charged lepton sector. Using (36), (37) and 

the relation I — I c 



Am 2 , n , 

t |°' , we find 



tan 6>i3 ~ 



Ami, tan 12 tan 23 



Using the current data [23], we find 



'13 



AmL l + tan 2 0i 2 



0.05-0.14 . 



(38) 



(39) 



For the central values 



Am 2 



sol 



7.9 • 10- 5 eV" 



Am 2 



atm 

2 , 



2.4 • 10" d eV 2 



sm" tfi2 = 0.314 , sin^ 23 = 0.44 , (40) 

(38) yields (0i 3 ) cent — 0.075. Now we can derive the (1, 3) element of the leptonic mixing matrix 
U l . A sizable contribution from the charged lepton sector comes only for i3 = \U l e3 \. Namely, the 
charged lepton (1,2) mixing angle Q epu ~ is relevant. Taking this into account, we obtain 







13 



\u. 



e3 I 



tan ^ 3 - , / — s 23 e 



m e 



AS-x) 



m. 



COS ^ 3 , 



(41) 



where angles and phases are defined in (29), (35)- (38). Since the phase 5 — x is unknown, we can 
calculate a range for \U l e3 \. Namely, from (41), 



for d\ 3 



for d\ 3 



0.05 , 0i3 = 0.01 - 0.09 , 

= 0.14, 0x3 = 0.08 - 0.2. (42) 
Thus, finally we get 6 

013 = 0.01-0.2. (43) 

The upper range in (43) is consistent with the current experimental bound [25], while future 
experiments [26] should be able to probe values close to 0.01. 



In summary we have proposed a relatively simple extension of minimal 5*0(10) by introducing 
two 'matter' 10-plets (and possibly one pair of vector-like matter 16 + 16 allowing us to predict the 
value of 0i3) plus some 5'O(10) singlets. Augmented by a i?-symmetry these states play an essential 
role in understanding the hierarchies in the charged fermion sector and in realizing bi-large mixing 
in the neutrino sector. The third leptonic mixing angle turns out to be naturally suppressed , in 
the range of 0.01 — 0.2. The /^-symmetry also implies 'matter' parity, so that the LSP is stable. 
It also forbids Planck scale suppressed baryon number violating dimension five operators. Since 
i?-symmetry has previously been shown to play an essential role in the construction of realistic 
S'O(IO) inflation models [27], [15], it would be interesting to try to merge the two approaches and 
also include leptogenesis. This will be attempted elsewhere. 

6 Numerically similar values have been recently predicted in [24]. 
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